HARDY-PETROVITCH-HUTCHINSON'S PROBLEM AND 
PARTIAL THETA FUNCTION 

VLADIMIR PETROV ROSTOV AND BORIS SHAPIRO 

To the memory of Vladimir Igorevich Arnold 

Abstract. In 1907 M. Petrovitch [15] initiated the study of a class of en- 
tire functions all whose finite sections are real-rooted polynomials. He was 
motivated by previous studies 112j of E. Laguerre on uniform limits of se- 
quences of real-rooted polynomials and an interesting result of G. H. Hardy [6]. 
(The term 'section' meaning truncation although quite inappropriate seems to 
be standard in the literature on this topic.) An explicit description of this 
class in terms of the coefficients of a series is impossible since it is deter- 
mined by an infinite number of determinantal inequalities one for each de- 
gree. However, interesting necessary or sufficient conditions can be formu- 
lated. In particular, J. I. Hutchinson [7] has shown that an entire function 
p(x) = ao + a\x + ... + a n x n + ... with strictly positive coefficients has the 
property that any its finite segment aix' + ai+ix' +1 + ... + ajxi has all real 

2 

roots if and only if for all i = 1, .... n — 1 one has 8 > 4. In the 

2 

present paper we give sharp lower bounds on the ratios — — 1 „ , i = 1,2, ... 
for the class considered by M. Petrovitch. In particular, we show that the 
limit of these minima when i — > oo equals the inverse of the maximal positive 
value of the parameter for which the classical partial thcta function belongs to 
the Laguerre-Polya class. We also explain the relation between Newton's and 
Hutchinson's inequalities and the logarithmic image of the set of all real-rooted 
polynomials with positive coefficients. 



1. Introduction 

In what follows we will use the terms 'real-rooted polynomial' and 'hyper- 
bolic polynomial' as synonyms. 'Hyperbolic polynomials' are inherited from par- 
tial differential equations. Consider the space V „ of polynomials of the form 
pix) = 1 + a\x + ... + a n x n with real coefficients. A polynomial p(x) € V n with 
all positive coefficients is called section-hyperbolic if for all i = 1, ...,n its section 
1 + a\X + ... + aix 1 of degree i is hyperbolic. Let A„ C V„ be the set of all 
section-hyperbolic polynomials of degree n and let A = (J A„ be the set of all 
section-hyperbolic polynomials. (Notice that by a result of E. Laguerre [H] a for- 
mal power series whose sections belong to A is an entire function lying in the 
Laguerre-Polya class.) We will call an entire function with positive coefficients all 
whose sections belong to A section-hyperbolic. The following question was posed 
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to the second author in April 2010 by Professors O. Katkova and A. Vishnyakova 
who attributed it to Professor I. V. Ostrovskii, comp. |14j . [2"T] . 

Problem 1 (Hardy-Hutchinson- Petrovitch-Ostrovskij^]) . For a given positive inte- 
ger i find or estimate 

f 

Denote by Pol n the space of all monic real polynomials of degree n and by 
£„ C Pol n the set of all such polynomials having all real and negative roots. 
Given a polynomial p(x) — ao + a\x + ... + a n x n , a n ^= define its reverted 
polynomial as P(x) = a^x n + a\x n ~ l + ... + a n = x n p(l/x). Obviously, reversion 
is a diffcomorphism between V n and Pol n . We say that a polynomial P(x) — 
x n + aix n + ... + a n is reverted section-hyperbolic if its sections x n + ... + a,iX n ~ l 
from the back are hyperbolic for all i — 1, ...,n. Reversion sends diffeomorphically 
the set A„ c Vn of all section-hyperbolic polynomials onto the set of all reverted 
section-hyperbolic polynomials. By a slight abuse of notation we denote the latter 
set by A„ C Pol n and will freely use both interpretations. (In fact, the second 
interpretation can be already found in the original paper j!5|.) Observe that the 

a 2 

quantities — t „ are preserved by reversion (up to a change of indices). 

Notice that since the natural projection 7r,; : A 1+ i — > A$ 'forgetting' the leading 
monomial is surjective one has 

2 2 

_ ■ f a i _ ■ f n 

m,i — «njpgA — m J?eA i+ i > 

i.e., to determine m,; it suffices to consider Aj+i. Moreover, by Hutchinson's the- 
orem, see [7], one has nij < 4. (Hutchinson's theorem was rediscovered 70 years 
later in [S].) Apparently Petrovitch knew that mj = 4, m 2 = ^ = 3.375 and 
m 3 w 3.264, see pp. 42-43 of [TS] but he made no explicit claim. On p. 331 of 
7J Hutchinson writes that the sequence {mi] should be strictly decreasing to some 
unknown limit of which it is known that it should exceed 2. 

The next result was proved by O. Katkova and A. Visnyakova around 2006. 

Theorem 1. For a given positive integer i one has 

rrii > 3. 

2 

Remark 1 . Observe that the quantities a °^ i can attain arbitrarily large values 
on A n so it is their minimal values which are important. Also notice that if one 
instead of A„ considers these quantities on the set £„ of all monic polynomials 
with negative roots, then the famous Newton's inequalities claim that for any given 
i = 1, ... , n — 1 one has 

(1) _^L_ > (n-» + l)(» + D > 



The order of the names is chronological in accordance with the years of their contributions 
to the topic under consideration. According to late Vladimir Igorevich Arnold (whose students 
we both were) if a mathematical concept is named after some mathematician this in most of the 
cases is due to the fact that this person has never considered this particular concept. 
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see e.g. [13] and comp. Proposition [6j Moreover, the equality in ([lj is attained 
exactly on polynomials with all real and coinciding roots. In particular, for ap- 
propriate choices of n and i the quantity a?/aj_iaj+i can be arbitrarily close to 
1. 

Below we solve the problem of Hardy-Petrovitch-Hutchinson-Ostrovskii by pre- 
senting 'explicitly' the entire function which simultaneously realizes all the above 
m,. Namely, consider the following sequence of polynomials pi(x) — 1 + x, P2{x) = 
1 + x + x 2 /A, P3(x) = 1 + x + x 2 /4 + x 3 /54, p 4 (x) = 1 + x + x 2 /4 + x 3 /54 + 
(69 — llv / 33)a; 4 /13824, ... given by the following inductive procedure p n {x) — 
p n -i(x) + A n x n , where A n is the maximal positive number such that P n {x) is 
hyperbolic. Denote by Poo{x) = 1 + x + Y^=2 A n x n . (This series appears on p. 42 
of |15j and one can show that p^ is an entire function.) 

The first result of this paper is as follows. 

Theorem 2. For any positive integer i one has 

a 2 A 2 



inf. 



peA - 



-lfli+l Aj_iAj + i 



where Ai are the above coefficients. In other words, m,i is attained at Pi+i € Aj+i 
or, equivalently, the function p^, (x) minimizes all simultaneously. Moreover, 
Pi+i is the unique (up to a scaling of the independent variable x) polynomial in 
Aj+i minimizing the quantity a|/aj_ia,+i. 

a 2 

By the above remark one can also conclude that i attains its minimum 

at the monic polynomial Pj+i € Poli+i which is the reverted polynomial to Pi+i- 
We call the inequalities of the form 

(2) a * > na 

O-i-lO-i+l 

with rrii defined above Petrovitch 's inequalities and the inequalities of the form 

a 2 

(3) > 4 

Hutchinson's inequalities, see footnote above. Using Theorem [T] one can show that 
rrii are algebraic numbers and calculate them on computer with an arbitrary pre- 
cision. Namely, the 10 decimal places of the first 17 mi's are as follows: mi = 
4, m 2 = f , m 3 = 2(69+ 8 1 1 lV55) « 3.2639552867, m 4 « 3.2403064116, m 5 w 
3.2351101647, m 6 « 3.2339623707, m 7 « 3.2337086596, m 8 w 3.2336525783, 
m g w 3.2336401824, mi « 3.2336374426, m u « 3.2336368370 m u w 
3.2336367032, m i3 « 3.2336366736, m u w 3.2336366671, m 15 « 3.2336366656, 
m 16 sa 3.2336366653, mi7 ~ 3.2336366652. Notice that for all i > 17 all have 
their first 10 decimal places coinciding with that of myj. Further calculations show 
that the next ten rrii have the same 10 decimal places as myj, they are monotone 
decreasing, and every second time the next decimal position stabilizes. Our next 
result confirms this behavior. 



Theorem 3. The sequence {mi},i = 1,2, ... is strictly monotone decreasing. 
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From Theorems [T] and [3] we get that the sequence {mj has a limit which we 
denote by . The main result of this paper is the explicit description of moo and 
the related entire function. To do this we define the formal power series ^(q, u) 
in the variables (q, u) which we by a small abuse of notation call the partial theta 
function 

00 

(4) = y2 q uj - 

3=0 

This function already appears on p. 330 of [7]. 

Remark 2. The standard partial theta function is usually defined by the series 

Q(q,u) = X)jlo(~l)"' < ?''^ uJ ' see e -§- ffl' 0' an d these two functions satisfy 
the obvious relation 

(5) $(q,u) = Q(q,-qu) 

which allows to translate their properties to one another. A number of beautiful 
identities which it satisfies was stated without proofs in Ramanujan's "lost" note- 
book which was found by G. E. Andrews after 50 years and who put significant 
effort in proving these identities. New results about the sum and product of partial 
theta functions can be found in e.g. [3]. It is also of interest in statistical physics 
and combinatorics, see [20] . 

In what follows we consider q as a parameter and u as the main variable. One 
can easily see that ^>(q, u) has a positive radius of convergence as a function of u 
if and only if \q\ < 1. If |g| = 1, then ^((^it) has a convergence radius equal to 1 
while for any such q with \q\ < 1 the function ^(q, u) is entire. Moreover for small 
positive q the series ^(q,u) considered as a function in u belongs to the Laguerre- 
Polya class C — V + , i.e. it has all real and negative roots, see e.g. [IT], Ch. 8. (The 
well-known characterization of these functions was obtained almost hundred years 

ago in [I7j.) Notice that for ^(q,u) the quotient a^/aj-iaj +1 
equals 1/q. 

Recall that Pi denotes the reverted polynomial for pi introduced above. (Since 
the constant term of pi equals 1 one gets that Pi is monic.) In Lemma [8] below 
we show that each Pi has all simple negative roots except for a single double root 
which has the minimal absolute value among all roots of Pi . Since each root of Pi 
is the inverse of the corresponding root of Pi one gets that Pi also has all real and 
negative roots except for a single double root which has the maximal absolute value 
among the roots of pi. Denote by Q the unique double root of Pi. For a positive 
integer i = 1,2,... define the scaled reverted polynomial Pi(x) = Pi(—(ix)/Pi(0). 
The scaling of is done in such a way that its double root is placed at —1 and its 
constant term equals 1. The main result of our paper is as follows. 

Theorem 4. The limit rrirx, = lim^oo m, exists and coincides with 1/q, where 
q > is the maximal positive number for which the series ^(q, u) belongs to C — P + 
as a function in u. Moreover, the sequence {Pi} of the scaled reverted polynomi- 
als converges to ^!(q,—ux), where u is the unique real double root of ^(q,u), see 
Figure [7} 

Notice that by ^ the function <d(q,u) belongs to C — V + exactly on the same 
interval of values of q as ^(q,u) does, namely, q e (0, q). The approximative value 
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FIGURE 1. <&(§» in the intervals [-40,0] and [-150,0]. (The 
value u of the negative double root of ^(q, u) with 10 decimal 
places is -7.5032559833.) 



of q with 10 decimal places is 0.3092493386, see Figure [T] We will later show that 
q is a root of transcendental equation (111. Notice that the constant l/q « 3.2336 
has earlier appeared in the papers and [10] , where the authors studied functions 
closely related to the above partial theta function. Namely, one of the main objects 
in paper [9] is the function 

OO 

9a{x) := V -p-, a > 1. 

3=0 

Theorem 4 of [9] claims that g a (x) has all hyperbolic sections, i.e. is section- 
hyperbolic if and only if a 2 > l/q. (Similar statements can be found in a recent 
preprint [5].) Again, notice that there exists a simple relation between g a {x) and 
0(g, u), namely 

and, therefore, the functions g^/q(u), ^>(q,u) and <d(q,u) belong to C — V + exactly 
on the same interval (0, q) of values of q. 

Furthermore, Theorem 2 of |10j claims the following. 

Theorem 5. Let f{x) = Y^jLo a i x3 ' i a j > be an entire function and S n (x) — 
Sj=o a j x3 t> e Us sections. Suppose that there exists a subsequence {nj}°° =1 C N such 

a 2 

that S nj (x) is hyperbolic for j — 1, 2, .... If there exists <5oo(/) = hm n _ i . 00 - — ^ — 
then for any positive integer m one has that X)j=o j!(vT~ )j 2 ~ * s hyperbolic. 

As a corollary of Theorem [5] one obtains that if fix) is section-hyperbolic (at 
least for all sufficiently large n) and 5oo(/) exists, then <5oo(/) > l/q- 

Finally let us mention a technical result of independent interest very much in the 
spirit of the modern study of amoebas of complex hupersurfaces and, in particular, 
of discrimininants. Denote by I/£ n C 1" (respectively LA n C K.") the images of 
E n C Pol n (respectively of A„ C Pol n ) under taking coefficientwise logarithms. 

Proposition 6. (i) The polyhedral cone given by Hutchinson's inequalities ^ (i.e. 
coinciding with the logarithmic image of the set of two-sided truncation hyperbolic 
polynomials) is the maximal polyhedral cone contained in LA n . The same cone is 
the maximal polyhedral cone contained in LE„ . 

(ii) The minimal polyhedral cone containing LA/~ is given by Petrovitch 's inequal- 
ities ([2]) while the minimal polyhedral cone containing LT, n is given by Newton's 
inequalities ([T]), see Remark [7J 
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Notice that Hutchinson's cone is, on the other hand, the recession cone of the 
logarithmic image of the set sign-invariant hyperbolic polynomials, see |16) . A fact 
similar to Proposition [6] is proven in Theorem F of |10) . 

Acknowledments. The authors want to thank Professors O. Katkova and A. Vish- 
nyakova of Kharkov National University for the formulation of the problem and 
their proof of Theorem 1 which they kindly allowed us to include in the present pa- 
per. They also gave us important hints and contributed to the proof of Lemma [T7| 
We are sincerely grateful to Professors A. Eremenko and A. Sokal and, especially, 
Professor G. E. Andrews for the valuable information about partial theta func- 
tions. The second author wants to acknowledge the hospitality of Laboratoire de 
Mathematiques, Universite de Nice during his visit in April-May 2011 when this 
project was carried out. 

2. Proofs 

We start with Theorem [T] 

Proof. Take some polynomial r n (x) — oq + a\x + ... + a n x n belonging to V n and set 
Tk(x) — a,Q + aix + ... + a k x k , k = 2, 3, . . . , n. By our assumption the polynomials 

2 

r k (x) are hyperbolic for all k. Set '■= ^—^-^k '■= — — — — >LzJ —. 

Let us fix an arbitrary k = 3,4, ... ,n. Denote by > x[ k ^ > x^ > ■■■ > x^ 
the zeros of r k (x). Using the Cauchy inequality with Ix^l 1 / 2 and ja^l 3 / 2 we get 

From the standard identities for (elementary) symmetric functions we get 

(fc) . (fe) . . (fc) flfc-i 
x\ + x\> + ... + x\> = = -7 fe , 

I (*0\2 , / (fc)\2 , | / (fc)\2 ( a k-l\ 2 a fc-2 2 o 

(x\ ') +(4 +... + {x\ ') = -2 =7|-27 fc 7 fc _ 1 , 

V a k J a k 

(x[ k ^M4 k) ) 3 +-.M4 k) ) 3 - - f^V +3 ^I^-3^ = -,1+Hlk-^lklk-^k-,. 

\ a k J a k a k a k 

Substituting these identities in the above inequality and dividing by j^.j k -i we 
obtain 

7^ - 47 fe _i + 37fc_ 2 > 0. 

Dividing the latter inequality by 7^-2 and using Sj's we get the following in- 
equality: 

(6) 44- 1 "44-i +3 >0. 

Since k = 3, 4, . . . , n is an arbitrary index, we get from (|6| the following system 
of inequalities 



(7) 



44-i -44-i +3 > 0, fc = 3,4,...,n. 
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Since ra(x) is hyperbolic we have 62 > 4. Suppose that the statement of the theorem 
is not true, and denote by j the smallest index such that Sj < 3, so that Sj-i > 3 
and Sj < 3 (j = 3, 4, . . . , n). We rewrite Q for fc = j in the form 

- A)8j_ 1 + 3 > 0. 

Since <5j ■ — 4 < 0, and Sj < <5j-i, the above inequality implies 

(« i -4)5 i + 3>0, 

whence <5j € (— 00, 1) (J(3, +00). By our assumption rj(x) is a hyperbolic polyno- 
mial, thus Sj 6 (—00, 1) is impossible. We conclude that Sj > 3. 

□ 

To prove Theorem [2] we need some preliminaries. Observe that rescaling of 
the independent variable x by an arbitrary positive constant acts on all spaces of 
polynomials we introduced above preserving the quantities of / a n -\a n . This action 
allows us to normalize <Xi = 1 in Pol n and analogously a\ = 1 in V n and, therefore, 
to reduce the number of parameters by one. Define V 1 as the space of polynomials 
of the form p(x) = 1 + x + a^x 1 + ... + a n x n and Pol n as the space of polynomials 
of the form P{x) = x 11 + x 71 ^ 1 + a 2 x n ~ 2 + ... + a n . Notice that reversion sends 
V n onto Pol^ an d that our main polynomials pi belong to V\ while their reverted 
polynomials Pi belong to Pol}. From now on instead of working in A„ C Pol n we 
will work in A* C Pol n which is the restriction of A„ to Pol}. The above group 
action carries our proofs from one space to the other. 

Define the standard embedding enrij n : Polj — > Pol}, j < n (respectively 
Polj — > Pol n ) given by multiplication of a monic polynomial of degree j < n 
by x n ~ 3 . Obviously, the image emj yn (Polj) C Pol} coincides with the coordinate 
subspace of all monic polynomials having all coefficients of degree less than n — j 
vanishing. Denote by T>j C Pol^j the standard discriminant consisting of all monic 
polynomials of degree j having at least one real root of multiplicity at least 2. Em- 
bedding T>j into Pol} using erriij let us define the discriminant T>j n C Pol n by 
taking the trivial (n — j)-dimensional cylinder over emij(T>j) along all coefficients 
of degree less than (n — j). Define A* C Pol^ and C Pol n as the restrictions of 

A„ and £„ to Pol}. Finally, consider the closure A* C Pol n of the set A* C Pol}. 

Lemma 7. (i) The set A n has a natural stratification of an (n — 1)- dimensional 

simplex with vertices at Pi = x n + x n ,P% — x" + x n + x"~ 2 /4, P3, P n . 

Different (n — 2) -dimensional (b oundaryj faces of A n belong to different Dj^ n , j — 

0, 1, 2, n, see Figure^ 

' — 1 1 1 

(ii) The natural projection Tr n by 'forgetting' the constant term sends A n onto A n _ 1 . 

(iii) Any polynomial in A n can be connected to P n by a smooth path along which 
all coefficients are non- decreasing. 

Remark 3. The original set A„ (respectively A„) is the cylinder over A* (re- 
spectively A„) obtained by the action of the group of rescaling of x by positive 
constants. 

Proof of Lemma [?| The first two statements are rather obvious and proved by in- 
duction. A n _ 1 is naturally embedded in the hypcrplanc a n = of Pol} using the 
multiplication of polynomials of degree n — 1 by x. Then A n is fibered over the im- 
age of A n-1 in Pol}^ along the constant term. To prove the third statement we show 




FIGURE 2. The domains £3 and A3 for the family x 3 + x 2 + ax + b. 
(Notice that £3 is the largest curvilinear triangle containing A3.) 

that the face T> n of the boundary of A n can be expressed as a n — a n (a2 1 a n -i), 
where (02, a„_i) € A n _ 1 . Moreover, for each i = 2, n— 1 one has da n /da,i > 
in the whole open A^_j. Indeed, denote by Xi < the roots of P(x). The double 
root is denoted by i n _i = x n . Since a, > for all i = 1, ...,n it will be conve- 
nient to consider them as elementary symmetric functions in the positive quantities 
—Xj. One has da n /da,i — Y^Zi(^ a n/9(—Xj))/(dai/d(— Xj)). The quantities a„, 
da n /d(—Xj) and dai/d(~Xj) are given by homogeneous polynomials with positive 
coefficients in all — Xj. This fact implies that the direction derivative of the func- 
tion a n is non- negative along any vector in A„_j with all non- negative coordinates. 
Using this statement together with induction on n we get that any polynomial in 
A* can be connected to P n by a smooth path with non-decreasing coordinates. □ 

Remark 4. In Lemma[8]below we will prove that each polynomial in A,\ has either 
simple negative roots or at most one double root (in which case it belongs to V n ) 
which is the rightmost among all roots of the considered polynomial. 

Proof of Theorem^ Given some reverted section- hyperbolic polynomial P(x) = 
x n + x n ~ 1 + a 2 x 2 + ... + a n of degree n > 4 consider the function k := a 2 l _ 1 /a n -2 a n- 
We want to show that m n —i — min^i k is attained at P n (x) which is the reverted 
polynomial to p n (x) defined in the Introduction. For fixed 0,2, a«-i the function 
k is minimal when a n is maximal in which case the polynomial P(x) belongs to 
T> n . Thus, we can restrict our consideration to A^ 3 P(x) € T> n . Since each P(x) 
can be connected to P n (x) by a smooth path along which each coefficient is non- 
decreasing it is enough to show that for all i = 2, n — 1 the partial derivative 
dn/dcLi is negative when k is restricted to T> n . There are three different cases to 
consider: (1) i < n — 2; (2) i = n — 2 and (3) i = n — 1. For i < n — 2 one has 

dn 8k da n a n-i da n 

dai da n dai an-20^ ddi 

Since all a.; > and da n /da,i > on V n by Lemma [7] one has that case (1) is 
settled. Analogously, we have 

Ok _ 2 a n + a n - 2 (dan/dan-2) 
Again since all > and da n /da n - 2 > on V n case (2) is settled. 
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Finally, one has 

Ok _ 2a n -i a^^dan/dan-i) a„_i(2a„ - a„_i(9a„/9a„_i)) 



, where 



9a n _x a„_ 2 a„ a„_ 2 a£ (a„_ 2 a2) 

ga« _ y> da n I dOn-l 
n_1 da n _i n_1 ^ 9(-a;,) / 

For z 7^ rt— 1 one has a„_i = fi + higi, where fi and ^ are homogeneous polynomials 
with positive coefficients depending on the positive variables —Xk, k ^ i. Therefore 
for i =^ n — 1 one has 

da n -i da n / da n - 2 da n 
(8) a„_i > —XiQi = —Xi— and a„_i — / — > —x 



d(-Xi) ' d(-Xi) I d(—Xi) d(-Xi) 

For i = n — 1 one has a„_i = —x n -\v + x\_ x w and 

d(-x n -i) 

Here v and w are given by homogeneous polynomials with positive coemcients in 
—Xj. (In what follows we will not need their explicit formulas.) Therefore, one has 

x n -i <9a„_i 



Thus 



2 9(-x„_i)' 



9a„ „ 1 \ -> 9a 



1 ^ 5a„ / n\ 



2a„ - a n-iQ^ < 2a « + 2 Xi d{-x.-) < a " I 2 - « ' - 

For the homogeneous polynomial a n of degree n we used Euler's identity 

ri— 1 ^ 

, d(-Xi) 

By the above argument any directional derivative dn/du is non-positive if u is 
an arbitrary vector in A*^ with all non-negative coordinates. Moreover since any 
polynomial P £ A* can be connected with P n by a smooth path with nondecreasing 
(and on some subintervals strictly increasing) coordinates we have that the value 
of K at P n is strictly smaller at any other such P ^ P n . In the bigger set A„ this 
means that only polynomials obtained from P n by scaling of the variable x can 
have the same value of K = m n _i as P n has. The result follows. □ 

Now we settle Theorem |3] 

Proof. We will use induction. The base of induction is that 4 = m\ > m 2 = ^p. 
Assume now that the statement is proved for TOj_i and we want to show that 
rrii-i > 771$. By Theorem [2] TOj_i is attained as the quotient af_i/ai-2d>i at the 
polynomial P$ which is a monic polynomial of degree i. Moreover up to scaling of 
x the polynomial Pi is unique in Poli where this minimum is attained. 

Set Pi(x) — Ylj=olj x ' > ■> li = !• The quotient 71/7270 coincides with mj_i. 
Given a polynomial R denote by the result of the fcth truncation of R from 
the back, i.e. the polynomial obtained by removing all terms of R of degree smaller 
than fc. 
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Consider a perturbation R(x) := Pi(x) + eQ(x), where e > and Q is a monic 
polynomial of degree i + 1. We choose Q such that for all k = 0, ... ,i — 2 the 
truncation has a root at the unique negative double root of . (Notice that 
Pi' 1-1 -' has a single negative real root which we do not have to worry about.) 

Setting Q := x n+1 + Y^j=o a o x ^ one can easnv see that the latter condition yields 
a triangular linear system (T) for the undetermined the coefficients a. a. Hence it 
has a unique real solution which we denote by Q* . 

For e > small enough all coefficients of the polynomial R*(x) — Pi(x)+eQ*(x) 
are positive. All roots of all polynomials R*^ k \ k = 0, .. .,i — 1, are real, distinct 
and negative. 

Moreover, for the perturbation R* the quantity 

(71 + eai) 2 /(72 + ea 2 )(7o + £"o) 
tends to m.j_i as e — > 0. Therefore mj < mj_i. To prove that the latter inequality 
is strict we argue as follows. 

The quantity 71/7270 does not change when one performs a linear change of 
the variable x. Perform such a change after which the polynomial R* (up to a 

constant factor) becomes x l+1 + x 1 H . The linear change and the subsequent 

multiplication by a positive number increase the coefficient of x i+1 and decrease 
the coefficients of x l for I < i. The latter tend to as s — > 0. 

Now consider Pj+i. The above means that one can find <5 > and a sequence 
of degree i + 1 reverted section hyperbolic polynomials {i?} remaining outside the 
ball 5,5 centered at Pj+i and of radius 5 and for which the quantity 71/7270 tends 
to rrii. Indeed, all coefficients of Pj+i are positive while some of the ones of {Z} 
tend to 0. 

One knows that the minimal value of the quantity 71/7270 in Poh+i is attained 
only at Pj+i up to a scaling. Therefore there exists 77 > such that for all reverted 
section hyperbolic polynomials from dBg their quantity 71/7270 exceeds m, + 77. 

On the other hand, similarly to what we did while proving part (iii) of Lemma [7] 
one can define a procedure of continuously changing a polynomial Z into the poly- 
nomial P,+i so that the quantity 71/7270 strictly decreases. The continuous de- 
formation intersects B$. Hence Bg contains degree i + 1 truncation hyperbolic 
polynomials whose quantity 71/7270 is at the same time bigger than rrii + rj and 
less than some number arbitrarily close to mj_i. Hence rrii < TUi—i. □ 

In Lemma [7] we proved that is a curvilinear (n — 1) -dimensional simplex 
with vertices Pi = a;™ + x™" 1 , P 2 = x" + x n_1 + x™~ 2 /4, P n . Next we describe 
a corollary of Theorem [3] and Lemma [7] about the behavior of the function k = 
a n-i/ a «-2 a n on n — 1 edges of this simplex connecting the most important new 
vertex P„ with already existing vertices Pi, P n —\. 

Corollary 1. Fori = 1, ...,n — 2 the restriction of k — a 2 l _i/a„_2a„ onto the edge 
ei. n of connecting Pi to P n is monotone decreasing from the value m„_j_i to 
m n _i. On the remaining edge e„_i n the function n decreases from +00 to m n _i. 

Notice that formally k is not defined at Pi, ...,P„_i so the claim that «(Pj) = 
nfin-i-l i n Corollary [T] should be understood as a limit. 

We now proceed with Theorem [4] whose proof requires a number of intermedi- 
ate steps. In fact, we prove a more general statement. Recall that the sequence 
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{pi}, i — 1, 2, ... was obtained starting from p\ = 1 + x by adding to the previous 
polynomial Pi—i the maximal possible term a^x 1 , > such that the resulting 
polynomial is still hyperbolic. This procedure can be equally well started from an 
arbitrary polynomial Si(x) of some positive degree d with all negative and sim- 
ple roots. Thus, we obtain the sequence {si}, i = 1,2,... with Sj(x) = s,_i(x) + 
5>d+i-i% about which we will prove that hmi_ >00 (ai) 2 /ai_iai + i = 1/q, see 



Theorems 13 and 14 If we consider the corresponding sequence {Si}, i = 1,2, ... 
of the reverted polynomials, then one can easily check that 

(9) S j (x)=x(S :i - 1 (x)-S j - 1 (Z j - 1 )), 

where £j-i is a point of a local minimum of Sj—i(x) at which it attains the largest 
value among all its local minima. (Notice that, in general, is not unique. 

However, the sequence {Sj} is well-defined since by definition the value 
is the same for all possible choices of £j'-i-) Formula Q is well-defined if degSj_i > 
2. In the exceptional case deg Si = 1 we set S2 = xS\. For j > 2 define Tj := Sj/x. 

Lemma 8. The following facts hold: 

a) For j > 2 the polynomials Tj have all negative roots. 

b) Exactly one of these roots (namely, the one at is a double root and the rest 
are simple. 

c) For j > 3 the point is the rightmost critical point of Sj-i; the critical values 
of Sj—i at all other local minima are smaller than Sj—i(£j-i), i.e. the absolute 
values of all other local minima are larger than the one at £j-i- 

Proof. Denote by > X\ > ■ ■ ■ > Xd the negative roots of Si and denote by t any of 
its critical points, where Si has a local maximum. (Here we assume deg Si > 3. If 
deg Si < 2, then for S2 and S3 the above claim can be checked directly and starting 
with 54 we can use the same argument as for the case deg Si > 3.) Suppose 
that Xi+i < t < Xi. Then the polynomial T^o;) = S2(x)/x — Si(x) — 
is hyperbolic, with all roots negative. Under our assumptions the root of ^2(2:) 
at £1 has multiplicity 2. In principle, the critical value Si(£i) might be attained 
at more than one minimum of Si, in which case the polynomial T2(x) has other 
double root(s) as well. W.l.o.g. we might assume that £1 is the rightmost of 
these local minima, where the critical value is maximal among all local minima. 
Then on the interval (£i,0), the polynomial S2 has a unique local minimum. We 
temporarily denote it by £ 2 and show that at this minimum the critical value is 
maximal among all minima of S2 thus justifying our notation. Indeed, if v = S^fe) 
one has |v| = max^^o) |x||Si(:z;)-Si(£i)| < |6.||Si(£i)|. As \t\ > |£i|, one obtains 

|s 2 (t)| = (|S'i(i)| + |s 1 (Ci)|)|t|>|SiKi)||ei| = H. 

On the other hand, there exist two roots /„+i < f v < of S2 such that t € 
(/i/+i) fv)- The polynomial S2 has a local minimum on (f v +i,fu) an d for the 
critical value r of S% at this minimum one has \r\ > IS^WI hence \r\ > \v\. 

Thus S2 can attain the largest value among all its minima only at £2- This 
implies that except one double root at £2, the polynomial Ts(x) — Ss/x has all 
its roots distinct and negative. The same argument proves the statement for all 
i>3. □ 

Remark 5. Lemma[8]implies that if the initial polynomial Si has the property that 
its rightmost minimum has the largest (i.e. having the smallest modulus) critical 
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value among all its minima, then the same property holds for the whole sequence 
{Sj} defined above. 

In what follows we will always use the latter assumption on Si . Lemmas [91 ITT] 



and 12 summarize further properties of the polynomials Sj constructed under this 



assumption. 

Lemma 9. For each positive integer k one has 6/£fc— l < 1/3- Moreover, equality 
takes place only for deg Si = 1 and k = 2. 

Proof. Assume that deg Si = d > 2 which implies that degS^ = k + d — 1. Set 
as above Tj,_i := Sk-i(x) — Sfc_i(6-i)- The quantity 6 satisfies the equality 
£fc^]fe-i(£fc) + ?fc-i(6) = 0. Substituting the variable x by — 6-i a; we ma y assume 
w.l.o.g. that 6-i = — 1. Hence 

[ ' 6 7W6) 6 + 1 £j h- a i ' 

where otj are the roots of 7fc_i smaller than —1 (listed in the increasing order). 
The equation —1/6 = 2/(6 + 1) has the unique solution 6 = —1/3. One can 



easily check that this is the solution of ( 10 ) only for deg S± — 1 and k = 2 in which 



case Si = (x + 2), 5*2 = x(x + 2) and S3 = a; (a; + l) 2 (up to rescaling). For larger 

S~ 4 i/(a-^ 



fc the presence of the additional summand X^=i 4 V(£fc — a j) m the right-hand 



side implies that the graphs of the l.h.s. and the r.h.s. of ( 10 1 are intersecting each 



other closer to the origin than —1/3. Indeed, the function —1/6 is increasing on 
(—1,0) while the functions 2/(6 + 1) and Y^j=i~ A V(6fe — a j) are decreasing there. 
Each of these functions takes positive values on (—1, 0). The functions —1/6 and 
2/(6 + 1) tend to +00 when their arguments tend to and —1 respectively. □ 

For a given initial polynomial Si as above define A rn = —S m -i(£ m -i), i.e. A m 
is the absolute value of the largest minimum of S m _i. 

Lemma 10. For I > m > 1 one has A t < A m (4|6 l -i|)'~"73 (i ~ m)(i ~ m+5)/2 . 

Proof. By Taylor's formula applied at 61-1 one has for x € (6u 0] that S m -i(x) — 
S m -i(£ m -i) = (x - £ m -i) 2 b(t x ), where t x <E (61-1,2;) and b = Sl^-i- Hence 
S m (x) = x(x - £ m -i) 2 b(t x ). 

The function b is non-decreasing on [£ m , 0] and for m > 2 it is strictly increasing. 
Indeed, b(t) is the second derivative of the hyperbolic polynomial <S m -i having all its 
roots smaller than 6i-i < £m- Therefore, 6(6-1) > for Z = 0, 1, . . . , deg S m -i — 
1. 

For m > 2 the quantity is an increasing function of x. Indeed, consider the 
functions F :— (x — £ m -i) 2 b(t x ) and G := (x — t; m -i) 2 b(t Xl ). For > x 2 > x\ > 
61-1 one has F(xi) = G{xi) and F(x 2 ) > G(x 2 ). Therefore b(t X2 ) > b(t Xl ) hence 

Consider the quantity |S m (:c)| = \x(x — £ m -i) b(t x )\. Its maximum on [61-1, 0] 
equals A m+1 = \S m (t; m )\. Set R(x) := \x(x - f m _i) 2 |. Hence 



max R{x) = R — — = — |£ m _ 
aeKm-i.o] V 3 / 27 

On the other hand, 



-1 
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A m+1 = \S m (U)\ = R(U)KkJ < R (^y 1 ) & ( f o) = 

4 4 

= ^|£m-x|£ m _i&(£o) = ^ICm-llAn- 

This is the required inequality for I = m + 1. To obtain it for / = m + 2 recall that 
\(m\ < ICm-i|/3, by Lemma|9] Hence 

4 i / 4 \ 2 

An+2 < 7^|£m|An+I - 3 ( 27 j ^m-lAn ■ 

Suppose that A/ < A m (A^ l m m )/^ l - m ^ l - m+5 ^ 2 . Then < |Cm-x|/3'" m and 

< < ^A m (|e m - 1 |/3^™)(4|e m - 1 |)'-'73 ( '- m)( '- m+5)/2 = 

= A m (4|e m -i|) i_m+1 /3 a " m+1)( '" m+6)/2 , 
which proves Lemma [T0| by induction on i. □ 

Lemma 11. For 1 < s < k — 2 one has sgn Sk(^ s ) — (— l) fe ~ s+1 . 

Proof. By definition S m +i(a;) = cc(S l m (x) + Hence for Z > m one gets 

; 



S,(a:) = a;'- m 5 m (a;) + £ ^ 



i-i+x 



As 5 m (Cm-i) = 0, one has Si(£ m _i) = E ,= m +i ^iC-i • The si g ns of the terms 



i-j+X 



in this sum alternate (because £ m _i < and A m > 0). By Lemma 10 their absolute 
values rapidly decrease and it is the sign of 4 m +i^ which defines the sign of 
Si(^ m —i). Indeed, compare this term with the quantity 



B := 



- E Ajl&n-l" l<^m+x|CL-"ll E 3(j-m-X)(j-m+4)/2 
j — m + 2 j— m+2 



The sum in the right-hand side is majorized by J2v=i ^ /3 = 4/23, so 

and sgn <Si(£ m _i) = sgn = ( — l)'~ m . Setting s = m — 1 we get the required 

statement. □ 

Lemma 12. One has £ k /£k-i > 0.2864887043. 

Proof. Use the notation in the proof of Lemma [9] On the interval (—1,0) the right- 
hand side of equation ( 10 1 is majorized by 2/(£ fc + 1) + Ylj=i ~ V(£fc + & ) . Indeed, 



one has aj £ (Cj'+Xj £j+2)- To majorize one can replace l/(£fc — ay) by l/(£& — £j+x) 
and then use part a) of Lemma [8j 

As > —1, one can further replace l/(£fc + &) by l/(3- 7 — 1). For j > 2 one 
has 3 J - 1 > 2(3 J '" 1 - 1) and further 
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^ (3J - 1) < ^ (33 ' - 1) < 2 + 8 + 26 + 80 + 121 < 16' 

j=i j=i 

where Y^jLs — 1) is majorized by twice its first term which equals 1/242. 

Therefore will be majorized by the solution of the equation 

1 2 11 

+ 1 + 16 

which belongs to (-1,0). The latter equals -0.2864887043... implying that 
> 0.2864887043. □ 

Remark 6. The above upper and lower bounds for £,k/£,k-i are quite close to 
one another and imply that the quantities decrease approximately as a falling 
geometric progression. 

Define $(x) := -1/x - 2/(x + 1), ijj{r) := Y^Lt r J '/(l - r J+1 ). The next result 
is central in the proof of Theorem [4] 

Theorem 13. (i) The limit A = linifc-^oo £k/£k—i exists; 
(ii) A is the unique solution of the equation 

(11) $(-A) = V(A), 

belonging to (0, 1). 

Proo/. Set Z := 0.2864887043... and r = 1/3. Lemmas |9| and |l2| imply that if A 
exists, then it belongs to Iq := [Zo>f*o]- As in the proof of Lemma [9] rescale the 
variable x to obtain £k-i = ~ 1- 

We construct a series of closed intervals /; := [k, ri\, where U < Zj+i < rj+i < 7*j, 
such that for each fixed i one has £k-i G f° r & sufficiently large. 



Consider equation (10) and set U := Ylj=i 4 l/(^fc ~ a j)- Recall that otj G 



(£j+i)£j+2) f° r j ' > 1 where aj are the roots of T^-i in the increasing order. We 
can decrease the value of U at every point of the interval (—1,0) by assuming that 
for these j one has aj — and then by requiring to be as small as possible. 
The last condition means that each ratio £ m /£ m - 1 equals Iq for all m = 1, . . . , fc — 2. 



For each fixed k the solution to (the modified as above) equation ( 10 ) will exceed 
— 1/3 because —1/3 was obtained in the absence of the sum U. These solutions 
increase with k because when k increases more and more terms are added to U. 
Denote by — r*i the limit of these solutions as k — > oo. We get r x < r = 1/3. 

It is clear that if the limit A exists, then it must belong to the interval [Zo>nJ- 
Moreover, all accumulation points of the sequence of solutions to the original equa- 



tion (10) when k — > oo belong to [?o, r i]- 

Analogously, the value of U increases on the whole interval (—1,0) if one sets 
a j = £7+2 U — 1) an d requires ^- +2 to be the maximal possible. In this case each 
ratio £ m /£ m -i equals n for all m = 1, . . . , k — 2. 



Hence for each fixed k the solution to (the modified as above) equation ( 10 ) 
will be smaller than —Iq because — Iq was obtained when these ratios were equal to 
ro > r±. Denote by — Zi the limit of these solutions as k — > oo. As above we get 
that l\ > Iq. Thus Z < h < r\ < r . The construction of all further quantities Z, 
and Ti follows the same pattern. 
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Let us show that when i — > oo the lengths of the intervals Ii tend to as fast as 
a falling geometric progression implying that their common intersection is a single 
point. In other words, there is only one accumulation point of the sequence of 
solutions to the initial equation (10) as k — > oo. 

To do this present (10) in the form <&(&) = U(x) and let k —> oo. The modified 
equations will have the form 



(12) $(x) = <p(ri,x) and $(x) = ip(k,x), 

where ip(r,x) := X^iVC^ + (l/ r )"')- The left (respectively right) equation in 
(12) has — (respectively —r i+1 ) as its solution on (0, 1). The series ip converges 
uniformly on [Zo,?*o] x [ — 1,0]. 

Each of the functions <p(ri,x) and (p(k,x) is decreasing on (0,1). For each 
fixed x £ (0,1) one has tp(ri,x) > ip{li,x). Therefore, the intersection points 
of the graph of $(x) with that of tp(ri,x) and ip(li,x) belong to the rectangle 
[-ri,-k] x [(p(k,-k),tp(ri,-ri)]. 

For x G [-r ,-l } one has l/x 2 > 9 and 2/(a; + l) 2 > ^ > 3, hence |$'(x)| > 12. 
Therefore 

\n+i - k+i\ < ^Ivin, -n) - (p(k,-k)\ ■ 

To simplify the notation we write r instead of rj and I instead of U. Set 



oo 



M := ip(r, -r) - tp(l, -I) = 2j(r J ' - F)/((l - r 3+1 )(l - . 

For each j there exists 0j £ (/,r) such that r- 7 — V = j0 J ~ (r — I) < jri~ x {r — I). 
As I < r < 1/3, one has (1 - ^ +1 )(1 - P +x ) > (2/3) 2 = 4/9. Thus < M < 
(9/^)jr j - 1 (r - I). Recall that ^^Li J> i_1 = 1/(1 - r) 2 < 9/4. Therefore, 

\n+i - k+i\ < \<p(n,-ri) - ^(Zi,-/»)| < 12 81 16 ( r » ~ ^) < Tl 2 ~ ■ 

This proves part (i) of the theorem. 

To settle part (ii) one has to observe that the solution — ij+i to equation (12) 
and the parameter — both tend to —A, while 



Therefore, A solves equation (11). □ 

Remark 7. The number A is defined as limit when fc — > oo of the critical points of 
the function xT}.{— £f.x) belonging to (0, 1). Equation (10) has a solution in every 
interval (aj,otj+i). For every fixed j — k there exists the limit as fc — > oo of the 
solution belonging to (ctj—k, ctj-k+i)- The proof is just the same as the one in the 
above theorem. Denote all these solutions by Ci > C2 > • • ■ ■ 

The quantities \aj \ are growing as a geometric progression. Therefore the infinite 
product 



W:= (ar+lJlJ^-^J 
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is an entire function of genus 0. Hence for every fixed j — k there exists ctj = 
lim/j^oo £Jj_fc. Define the entire function V by the relation 

(13) V(x) := J W[t)dt j J W(t)dt 

Then the above quantities 2j are the zeros of the function V. Notice that with the 
above normalization one gets V(0) = 1. 

Theorem 14. The quantity A^ n / (A m -iA m+ i) tends to 1/A as m — > oo. 

Proof. By Taylor's formula of order 2 applied at £ TO _i the polynomial T m _i has the 
form T(x) = (x-£ m _ 1 ) 2 S^_ 1 (£ m _ 1 + 0(x-£ m _iX:E-&n-l))> where 0(x-£ m -i) € 
(0, 1). Notice that for x > £ m _i, 0(x — £ m -i) is uniquely defined. Indeed, S'J n _ 1 is 
increasing for x > £ m — l as the second derivative of a hyperbolic polynomial with 
all roots smaller than £ m _i. 

Thus A m _! = T m _i(0) - &_iS£_i(fm-i + 0(-£m-i)(-£m-i)). Notice that 
A rn = max [Sni _ 1)0 ] |x(a;-C m _i) 2 S'^_ 1 (^„ l _i + - £ m _i)(:r - £ m _ x ))|. The maxi- 
mum is attained at £ m and £ m /£m-i — > A as m — > oo. Thus 

^-m , (Cm — £m-l) 2 <Sm-l(£m-l + ^(£m — £rra-l)(£m — £m-l)) 



e 



The first quotient to the right (denoted by Y m ) tends to (A — l) 2 . Denote by F m the 
second quotient. In the same way one shows that A m+ \/A m = £, m +iY m +iFm+i- 
Hence 

A m £,m Y m F m 



A m —lA m +i £m+l Fm+1 

with ]im m _ ) . 0o (£ m /f m+ i) = 1/A and lim I „_ >00 (y;„/y; n+ i) = (A - 1) 2 /(A - l) 2 = 1. 
So to prove the theorem we need to show that lim m _ i . 00 (i 7 ' m /i 7 ' m+ i) = 1. 

The latter statement follows from the above remark. Indeed, set x >-> \£ m -i\x 
(resp. x i y |£ m |x). Then the function T m —i (resp. T m ), in the limit as m — > oo, 
becomes the function V defined by ( |13[ ). The numerators of the quotients F m 
and F m+1 tend to V"(— 1 + 8(— A + 1)(— A + 1)) while the denominators tend to 
V"(-l + 9(1)). Hence F m /F m+1 -> 1 as m -> oo. □ 



Proposition 15. 27ie function V defined by (13) enjoys the following properties: 

(i) V(0) = 1, V(-l) = V'(-l) = 0; 

(ii) V belongs to the Laguerre-Polya class £ — V + ; 

(hi) V satisfies the functional relation: V(x) = 1 + xV(\x)/V(— A). 

Notice that the latter relation implies that for any choice of A one has V(0) = 1 
and V(— 1) = 0. On the other hand, V'(— 1) = is an additional condition which 
together with (ii) determines A. 

Proof of Proposition \15\ Part (i) follows from the definition of the function V, see 
Remark [7j To prove part (ii) notice that the functions V and W are limits of 
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sequences of hyperbolic polynomials Hk with negative roots. If the roots are num- 
bered in the order of increasing absolute values, then for every fixed j the root ctj 
has a finite limit when k — ¥ oo. 

The modules of the roots increase faster than a geometric progression with ratio 
2.6. Hence the sequence of polynomials is uniformly convergent on any compact 
set O. Indeed, consider the product Ili^ArU — x/ai). The module of its logarithm 
is majorized by C^2°^ N \x/cti\ (where C > depends only on the set Q) which 
is arbitrarily and uniformly on small if one chooses N sufficiently large. To see 
this notice that | ln(l + y)\ < \y\ + \y\ 2 + \y\ 3 + • • • , where y — —x/ai, for i large 
enough the sum (which equals |y|/(l — \y\)) is smaller than 2\y\. The limit of such 
a sequence belongs to the class C — V + by definition. 

To prove part (iii) one has to recall the equality S m+ i(x) = x(S m (x) — S m (£ m )) 
or, equivalently, T m+ i(x) = xT m (x) — S, m T m (^ m )- The function V is the limit when 
m — > oo of the polynomials T m (— £ m —ix) multiplied by a nonzero constant so that 
T m (0) — 1- Rescaling sends the double root of T m (x) to —1. □ 

In order to finish the proof of Theorem [2] we study the problem which entire 
functions satisfy the properties given in Proposition |15| The following definition is 
crucial for our further considerations. Recall that <f(q, u) = Y^jLo Q^'* We say 
that a pair (q, u) is critical for ^(q, u) if \q\ < 1 and ^(q, u) as a function of u has 
a double root at u. 

Theorem 16. There exists an analytic in a disk \x\ < r, r > 1 function V(x) 
satisfying the relation 

(14) V(x) = -uq(xV(qx)+V(-q)) 

for some m6C* and \q\ < 1 as well as the boundary condition V'(—l) = if and 
only if the pair (q, u) is critical. 



Notice that (14 1 is exactly the relation (iii) of Proposition 15 with an undefined 
scalar factor u. 

Proof of Theorem \l()[ Assume that a function V(x) analytic in some disk \x\ < 
r, r > 1 satisfies the relation V(x) = —uq(xV(qx) + V(—q)) for some fixed ti ^ 
and q ^ 0. W.l.o.g. we can assume V^(0) = 1 which is equivalent to (3q = 1. 
Substituting the power series V(x) — fii^ m ^ ne latter relation one gets 

oo oo 

"^^Pjxi = — quV{— q) — qux^^ f3jq 3 ' x^ . 

j=0 0=0 

Comparing the coefficients at equal powers in the latter relation we get the system 
of equalities 

(15) -quV(-q) = 1, and 



Pi = -qu, f3 2 = -q U/9i, k = -q ul3 2 



, [3 k = -q k u(3 k - 



implying that flj — q\ 2 >{— u) 3 , j = 1,2,3,.... Substituting these coefficients in 

V(x) we get V(x) = X)j!lo9^ 2 ( — ux ) 3 ■ In terms of the partial theta function 
$(g, u) given by Q one gets V(x) = ^(q, —ux) and condition (15) takes the form 
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ty(q,qu) = — ^j. Let us show that it is equivalent to \t(g, u) = 0. Indeed, expand- 
ing V(q,qu) = -± we get qu^oQ^ W = 1+Ej=o9^ W +1 = 
0^1 + E*=o9^ 

u i+i = ^(q,u) = 0. Notice that for the power series 
expressing V[x) to have a convergence radius exceeding 1 it is necessary and suffi- 
cient to have \q\ < 1 in which case V(x) is entire. Now we use the last boundary 
condition V'(-l) = 0. With V(x) = ^(q,-ux) we get V'(x) = ~ud^{q,-ux), 
where d stands for the partial derivative w.r.t. second argument. Finally, 
V'{— 1) = — u^' u (q, u). Thus one gets the system 

[*(«,«) =0 f*(g,«) = 

\-tt« / t ,(g,u) = \<(g,u)=0 ' 

since u = is never a solution of ^(q, u) = 0. 

An arbitrary solution (q, u) of the latter system is exactly a critical pair in the 
above definition, i.e. q is such that the function ty(q,u) as a function of u has a 
double root at u. □ 

Denote by g g (x) = Y^k=o1 k2x > < q < 1 the partial theta-function (in this 
form), and by S n (q,x) = J2k=o1 k xk ^ s ntn Taylor section. We will need the 
following lemma. 

Lemma 17. For every real q € (0, 1) there exists a number to £ N such that for 
all n > 2m + 2 i/ie number of non-real zeros of S n (q, x) is not greater than 2m + 2. 

Proof. We will use the following well-known identity for \q\ < 1 

l_ q 2k 



fc=l y fe=l 
(see e.g. [TB], Chapter 1, Problem 56). By this identity we have 



1 + 2^(-l) fc g fe " > 0, 5 6(0,1). 

k—l 

Thus for every q € (0, 1) there exists m = 2s + 1 € N such that the inequality holds 

m 

(16) l + 2j](-l)V 2 > 0. 

fc=i 

Then for every n > 2m + 2 we have for all /c, to + 1 < k < n — m — 1: 



1 n n 

(-l) fc S n (g, -^) = £(-l)''-y a - aw = <T fe2 v -' 



j=Q j=0 

(k—m—1 k-\-m n 

(-iy- k q u - k)2 + £ Mr¥ J ~ fe)2 + E (-iy- k q u - k)2 
j—Q j—k—m j — fc+m+l 

= : Si + S2 + S3. 

We see that the summands in Ei are alternating in signs and their moduli are 
increasing. So the sign of Si coincides with the sign of the {k — m—1)- th summand. 
Therefore, signSx = (-i)k-m-i-k = ^_^- m _i = j-_ 1 ^_ 2s -2 = L Analogously 

the summands in S 3 are alternating in signs and their moduli are decreasing. So 
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the sign of S3 coincides with the sign of the (k + m + l)-th summand. In other 
words, sign S 3 = (_l)k+m+i-k = _ (_!)2 S +2 = 1 Therefore Si > and 



S 3 > 0. By (|16j) we have 

fc+m / m \ 

S 2 = q- k2 {-iy- k q (j - kf = q~* 1 + 2^(-l) fe g fc2 > 0. 

j — k—m \ k—1 / 

Therefore for every k, m + 1 < k < n — m—1, we have (—l)S n (q, — -J^ ) > 0. Thus 
for n > 2m + 2 the polynomial S n (q, x) has not less than n — 2m — 2 real zeros (and 
the number of non-real zeros of S n (q, x) is not greater than 2m + 2). □ 

Corollary 2. For every real q £ (0,1) the functions g q {x) and ^{q,x) have a 
finite number of non-real zeros. Moreover, the number of non-real zeros is a non- 
decreasing function of q. 

Proof. Since ^(q,x) is obtained from g q {x) by rescaling of x it suffices to consider 
g q {x) only. To prove the first statement fix an arbitrary q € (0, 1). By Lemma [l7| 
we know that there exists m such that the number of non-real roots of any Taylor 
section S n (q, x) for sufficiently large n does not exceed m. Assume that the function 
g q {x) has I > m non-real roots. Take small circles surrounding these roots and not 
intersecting the real axis. By the Hurwitz theorem all Taylor sections with large 
n should have exactly / roots in the union of disks bounded by these I circles. 
Contradiction. 

To prove the second statement consider the sequence of Taylor sections. We 
prove that the number of real roots of any section S n (q,x) and g q (x) itself is a 
monotone non-increasing function of q € (0,1). Indeed, the sequence q n is a 
complex zero decreasing sequence (CZDS) for any q £ (0, 1), see e.g. [4j. Thus for 
any < q\ < q2 < 1 one has that to obtain the section with the value of parameter 
qi from that of qi one has to multiply the coefficients of the former by (91/92)™ • 
Since the latter sequence is CZDS the result for sections follows. For g q {x) the 
same argument applies since it has only finitely many non-real zeros and is of genus 
0. □ 



Finally to finish the proof of Theorem [2] notice that by Proposition [15] and 
Theorem 16 the function V(x) satisfies the functional relation (iii), V'(— 1) = and 
belongs to C — V + which implies that A = q and the function V equals ^{q, —ux) 
where u is the double root of the function ^(q, x). □ 



Remark 8. A somewhat mysterious equation (111 describes the set of all criti- 
cal points of the limiting function ^(9, — ux) of which u is the only critical point 
belonging to the interval (0, 1). 



Notice that ^(q, x) belongs to C — V + if and only if q 6 (0, q\ where q is the 
constant appearing in Theorem|4j comp. Theorem 4 of [5]. The following additional 
statements are very plausible but since we do not need them we have not put enough 
effort in proving them. 

Conjecture 1. (i) Enumerating the negative roots of^(q,u) as —Ti(q) < —r 2 (q) < 
—rs(q) < ... in the order of increasing absolute value we getliuiff^^ r " +1 ( ^ = l/q. 
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(ii) Analogously, if —r[ < — r' 2 < —r' 3 < ... stands for the real negative roots 
of ^f' u (q,u) or, in other words, for the real critical points of ^(q, u) we get 

limjv^ = l/ q . 

(iii) Finally, let ci,C2,..., be the sequence of the critical values of^(q,u) attained 

2 

at the real critical points r^,r 2 , then limjv->.oo CJV " CN x = 9- 

We call the set of all < 1 such that the entire function ^/(q,u) has a 

double root (denoted by u{), i.e. the pair (qi,Ui) is critical spectrum S of partial 
theta function ^(q, u). We propose the following conjecture about the properties 
of SSf(q,u) for q G (0, 1). 

Conjecture 2. The spectrum & contains infinitely many positive values < q — 
q\ < q_2 < ■■■ < qN < ■■■ < 1, i-e. there exists an infinite sequence {qi}, i = 1,2, ... 
of numbers in (0, 1) such that for each positive integer i the function ^(g,, u) has a 
negative double root in the variable u. 

Some numerical evidence for the validity of Conjecture [2] was newly obtained by 
two talented high school students A. Broms and I. Nilsson. Namely, they calculated 
the first 25 values of qi with 12 decimal places. Their list with 6 decimals is as fol- 
lows: 0.309249, 0.516959, 0.630628, 0.701265, 0.749269, 0.783984, 0.810251, 0.830816, 
0.847353, 0.860942, 0.872305, 0.881949, 0.890237, 0.897435, 0.903747, 0.909325, 0.914291, 
0.918741, 0.922751, 0.926384, 0.929689, 0.932711, 0.935482, 0.938035, 0.940393, see 
Fig. [| 

Problem 2. Assuming the Conjecture [2] holds what is the asymptotics of the 
sequence {qi} when i — > oo ? 

Let us also mention the following question posed by Professor A. Sokal. 

Problem 3. Is it true that & is empty within the open disk \q\ < q ? 

We finally prove the remaining Proposition [6j Denote by Pol+ C Pol n the set 
of all monic degree n polynomials with all positive coefficients. It contains S n , the 
set of degree n polynomials with all roots real negative. Let p{x) = a n x n H — • + clq 
be a polynomial. 

Lemma 18. For any k — 1 , 2, . . . , n — 1 there exists a polynomial p G Pol^ , p 
for which one has af > 4a^iai + i for i ^ k and a\ < Aak-iOk+i- 

Proof. Fix the triple of coefficients (<ife-i, dfe, cik+i) such that a\ < Aak-ia-k+i- 
Hence the polynomial g := x k ~ 1 (ak+ix 2 + a^x + ak-i) has a (k — l)-fold root at 
and a complex conjugate couple. 

For i > k + 1 (resp. for i < k — 1) set a% — biS l ~ k ~ 1 (resp. = bie k ~ 1 ~ t ), where 
£ G (0, 1]. Choose the coefficients bi > such that the inequalities af > 4ai_ia i+ i 
hold for i ^ k and e = 1. One can do this consecutively. E.g. one chooses bk+2 > 
sufficiently small so that a\ +l > 4afc&fc+2, then bu+s such that b\ +2 > Abk+30,k+i 
etc. Then in the same way etc. 

If the inequalities af > 4a.i_iai + i (i ^ k) hold for e = 1, then they hold for any 
£ G (0, 1] (to be checked directly). 

For £ small enough the polynomial p is a perturbation of the polynomial g. Fix 
two circles centered at the complex roots of g and not intersecting the real axis. For 
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Figure 3. The first 25 values of 



e > small enough one has \p — g\ < \g\ on these circles. By the Hurwitz theorem 
p has a root inside each of them. Hence p $ S n . □ 

Proof of Proposition^ To prove (i) notice that by Lemma 5 of [16] Hutchinson's 
cone is the minimal polyhedral cone containing the set of the so-called sign- 
independently hyperbolic polynomials and is, on the other hand, contained in 
A„ C S„. (A sign- independently hyperbolic polynomial is a hyperbolic polyno- 
mial with all positive coefficients and such that any sign change of its coefficient 



results into a hyperbolic polynomial, see |TB] .) By Lemma 18 we see that an arbi- 
trarily small parallel translation 'outward' of any of the hyperplanes defining the 
logarithmic image of Hutchinson's cone results into getting outside the logarithmic 
image of the largest set £„. Therefore, the logarithmic image of Hutchinson's cone 
is the largest polyhedral cone contained in LA n and, analogously, in LY> n . 

To prove (ii) observe that Theorem [2] can be interpreted as follows. Consider 
Hutchinson's cone in the space Pol 1 ^ and its image under the logarithmic map. Then 
it has a unique apex, i.e. the vertex where all inequalities become equalities. Look 
for a parallel translation of the logarithmic image of Hutchinson's cone containing 
the whole LA n . Then if you take the parallel translation when this apex is placed 
at the logarithmic image of P n , then the whole LA n is covered. Since the translated 
cone and LA n still have a common point this position is minimal for containment. 
Exactly the same argument using the known properties of Newton's inequalities 
tells us that placing the apex at the logarithmic image of (x + l/n) n does the 
job. □ 

3. An interesting iteration scheme 

The main ingredient in the proof of Theorem [4] is construction of the sequence 
{Sj} starting from a hyperbolic polynomial Si with all simple roots the rightmost 
of which is at the origin. The next polynomial is obtained from the previous one 
by subtracting its least in absolute value minimum followed by multiplication by 
x. We have shown that after appropriate scaling the limiting entire function is a 
specialization of a partial theta function and also that the critical point at which 
the minimum is located asymptotically stabilizes. These results naturally lead to 
the question about what happens if we consider a similar iterative procedure where 
the point at which we take a value to subtract is fixed from the beginning. A more 
detailed consideration leads to the following natural set-up. 

Given an initial analytic function f\{x) defined (at least) in a small open neigh- 
borhood of the interval [—1,0] on the real line and a number < q < 1 define the 
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sequence {fj} given by 

Obviously, fj(x) will be well-defined and analytic in the same neighborhood of 
[—1,0] unless — q is a root of fj-i(x). For generic choices of /i the latter circum- 
stance never happens. One can easily check that for a positive integer j if fj is 
well-defined, then it satisfies the normalization conditions /,(0) = 1 and fj(—l) = 0. 
Looking at the first part of the proof of Theorem [l6| before we use the additional 



condition V'{— 1) = one can see that the fixed points of (17), i.e. the analytic 
functions satisfying on [0, 1] the functional relation 

F(x)=(l ■ xF{qx) ^ 



F(-q) J ' 

with q G (0, 1) fixed are exactly of the form F(x) — ^(q, —ux) where u is one of 
the real roots of the equation ^>(q, u) = 0. Here u is the variable. 

Thus for any fixed q € (0, 1) the iteration scheme (JTYJ) considered as a self- map 
of an appropriate space of analytic in an neighborhood of [—1,0] functions has 
countably many fixed points. At the moment it is by no means clear what local 
properties these fixed points have. For example, for which q and which of the above 
fixed points are repelling/attracting? Under which additional assumptions on f\ 
the sequence {fj} obtained via scheme p7| converges? 

Our computer experiments and Theorem [4] suggest that the following statement 
should be true. 

Conjecture 3. If0<q<q and the initial function of the form fi = (x + l)Q(x) 
where Q{x) is a hyperbolic polynomial with all negative roots smaller than —1, then 
the polynomial sequence {fj{x)} converges uniformly on [—1,0] with any number of 
derivatives to the function ^(q, —u{q)x), where u{q) is the negative solution of the 
equation ^>(q,u) = with the minimal absolute value. 

Remark 9. Apparently under the condition < q < q the attraction domain 
of the latter fixed point is much larger than /i of the form given in the above 
conjecture. In particular, iterations started with f\ — sin7ra: converge very quickly 
to the same limit. On the other hand, if q > q numeric experiments show that 
iterations typically diverge, see Fig 3. It might be that for q > q all the fixed points 



of (171 become repelling. There are superficial similarities of the scheme (171 and 
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Figure 5. Convergence of iterations for q = 1/4 (left) and the 
critical q (right). 




Figure 6. Divergence of iterations for q = 1/2. The number of 
iterations on the left is 10 and on the right is 50. 



the famous logistic map in dynamical systems which also depends crucially on the 
value of the additional parameter q. 
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